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ABSTRACT 

We present a comprehensive study of the p Ophiuchi molecular cloud that addresses aspects of the 
physical structure and condition of the molecular cloud and its photodissociation region (PDR) by 
combining far-infrared and submillimeter-wave observations with a wide range of angular scale and 
resolution. We present 1600 square arcminute maps (2.3 sq. pc) with 0.1 pc resolution in submillimeter 
CO(4 — > 3) and [C I]( 3 -Pi — > 3 Po) line emission from the Antarctic Submillimeter Telescope and 
Remote Observatory (AST/RO), and pointed observations in the CO (7 -> 6) and [C I]( 3 J P 2 -> 3 -Pi) 
lines. Within the large scale maps, smaller spectral line maps of 3000 A.U. resolution over ^90 square 
arcminutes (0.2 sq. pc) of the cloud in CO, CS, HCO + and their rare isotopomers are made at the 
Heinrich Hertz Submillimeter Telescope Observatory (HHT) in Arizona. Comparison of CO, HCO + , 
and [C I] maps with far-infrared observations of atomic and ionic species from the Infrared Space 
Observatory (ISO), far-infrared and submillimeter continuum emission, and near-infrared H2 emission 
allows clearer determination of the physical and chemical structure of the p Oph photodissociation 
region (PDR), since each species probes a different physical region of the cloud structure. Although a 
homogeneous plane parallel PDR model can reproduce many of the observations described here, the 
excitation conditions needed to produce the observed HCO + and [O I] emission imply inhomogeneous 
structure. Strong chemical gradients are observed in HCO + and CS; the former is ascribed to a local 
enhancement in the H2 ionization rate, the latter is principally due to shocks. Under the assumption of 
a simple two component gas model for the cloud, we find that [C II] and [C I] emission predominantly 
arises from the lower density envelopes (10 3 to 10 4 cm -3 ) that surround denser cloud condensations, or 
"clumps". The distribution of [C I] is very similar to C 18 0, and generally consistent with illumination 
from the "far" side of the cloud. A notable exception is found at the the western edge of the cloud, 
where UV photons create a PDR viewed "edge-on" . The abundance of atomic carbon is accurately 
modeled using a radiation field that decreases with increasing projected distance from the exciting 
star HD 147889 and a total gas column density that follows that of C 18 0, decreasing toward the edges 
of the cloud. In contrast to conclusions of other studies, we find that no non-equilibrium chemistry is 
needed to enhance the [C I] abundance. Each spectral line is traced to a particular physical component 
of the cloud & PDR. Though CO rotational line emission originates from both dense condensations 
and diffuse envelopes, the millimeter-wave transitions mostly find their origins in envelope material, 
whereas the high-J submillimeter lines stem more from the dense clumps. Submillimeter HCO + and 
infrared [O I] and [C II] emission indicate clump surface temperatures of 50-200K, an ultraviolet 
radiation field with Itjv ~ 10 — 90, densities of 10 5 to 10 6 cm~ 3 , and interior temperatures of < 20K. 
This study highlights the value of large-scale infrared and submillimeter mapping for the interpretation 
of molecular cloud physical and chemical structure, and important future observations are highlighted. 
Subject headings: ISM: clouds — ISM: molecules — ISM: atoms — ISM: individual(Rho Oph) — radio 
lines: ISM 



1. INTRODUCTION 

With the advent of millimeter-wave spectroscopy, 
large-scale surveys of the 115 GHz CO ( J = 1 — > 0) line 
(e.g. Dame et al. 2001) have revolutionized our under- 
standing of molecular clouds, the birthplaces of all stars. 
Yet, in terms of clouds' physical structure and evolution, 
many elementary questions remain shrouded in uncer- 
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tainty. Illuminating these fundamental issues will enable 
a more complete understanding of how stars form (e.g. 
Bally et al. 1987) in a Galactic context. Sadly, numer- 
ous observational challenges frustrate these efforts. For 
example, optically-thin spectral line diagnostics of CO 
that probe the contents of the entire cloud entangle all 
structural cloud components along a given line of sight. 
One solution to this problem is to constrain one's view 
lM<5)firaii3<iiMe*iuctural regions of a cloud; for example a 
cloud's "surface" , where ambient ultraviolet (UV) pho- 
tons have not been entirely absorbed by dust grains and 
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(e.g. H2) molecules, and therefore illuminate (and dom- 
inate) physical structure and chemical processes. 

Photodissociation regions (PDR's) are these boundary 
interfaces, representing the transition between atomic 
and molecular phases of the interstellar medium (ISM). 
Most models invoked to explain the physics in these re- 
gions assume molecular clouds of uniform properties sur- 
rounding the UV radiation source, typically a neighbor- 
ing early type star (Ticlcns and Hollcnbach 1985; van 
Dishoeck & Black 1988; Burton, Hollenbach, & Tielens 
1990; Hollenbach et al. 1991; Le Bourlot et al. 1993; 
Koster et al. 1994; Kaufman et al. 1999). UV pho- 
tons with energies greater than the hydrogen ionization 
threshold of 13.6 eV create an H II region surrounding 
the energetic star(s). UV photons with energy less than 
this threshold escape from the H II region into an adja- 
cent molecular cloud, where photodissociation of domi- 
nant molecules like H2 and CO takes place, as well as 
photoionization of atoms with an ionization threshold 
lower than that of hydrogen - such as carbon, which 
is abundant and ionizes with photon energies as low as 
11.26 eV. The combined effect of CO photodissociation 
and C° photoionization by UV photons produces a lay- 
ered structure of C + /C°/CO surrounding the H II region 
and extending into the molecular cloud. This basic pic- 
ture of the C + /C°/CO interface structure, as modeled for 
homogeneous cloud structures, is not in very good agree- 
ment with observations of edge-on PDRs (i.e. M17); in 
general it is unable to account for the extended distri- 
bution of [C II] and [C I] emission seen in these regions 
(Stutzki et al. 1988, Keene et al. 1985). However, models 
of inhomogeneous photodissociation regions have proven 
capable of accounting for the extended [C II] and [C I] 
distribution observed in some PDR's (Stutzki et al. 1988; 
Boisse 1990; Tauber & Goldsmith 1990; Howe et al. 1991; 
Pineau Des Forets, Roueff, & Flower 1992; Meixner and 
Tielens 1993 [MT93]; Spaans 1996; Plume et al. 1999). 
The inhomogeneity, or "dumpiness" , of molecular clouds 
is well recognized (e.g. Blitz & Williams 1987); it seems 
only natural that such clouds' UV-illuminated surfaces 
should bear similar physical traits. 

An ideal observational target for such studies is the p 
Ophiuchi molecular cloud. It is one of the nearest (« 130 
pc) and most active star-forming clouds and has a rich 
and varied environment that provides the observational 
basis for the current "standard model" of low-mass star 
formation (Shu 1977; Shu, Adams, & Lizano 1987). The 
region hosts three hot B-type stars, although HD 147889 
appears to dominate the heating and excitation of the 
cloud's atomic and molecular "surface" (Liseau et al. 
1999). HD 147889 is the most-embedded (E B -v=l-09; 
Ay ~ 4) optically-visible star that has been studied by 
high resolution ultraviolet and optical spectroscopy; the 
line of sight passes through the diffuse atomic and molec- 
ular "halo" of the p Oph cloud (cf. van Dishoeck & de 
Zeeuw 1984; Crutcher & Chu 1985; Cardelli & Waller- 
stein 1986). Because of its deep obscuration and periph- 
ery to the main molecular cloud, HD 147889 is thought to 
lie just beyond the molecular cloud, illuminating it from 
behind. Modeling of the 158 pm [C II] line intensity as a 
function of projected distance from HD 147889 led Liseau 
et al. (1999) to the conclusion that the UV-illuminated 
PDR back surface of the p Oph cloud is approximately 
concave. 



Study of such a varied and complicated environment 
as p Oph is challenging. High spectral and angular res- 
olution is required to disentangle the complicated struc- 
ture of actively star-forming environments, but simul- 
taneously the overall structure of the cloud only be- 
comes apparent with studies of large angular extent. 
We have attempted to bridge these two disparate ex- 
tremes by exploiting the powers of far-infrared and sub- 
millimeter spectroscopy on both large and small tele- 
scopes. In § 2 and § 3, we present and interpret large- 
scale (40' x 40') maps of the molecular cloud in submil- 
limeter CO(J = 4 -> 3) and [C I] 3 Pi -> 3 P emission, 
with 25,000 A.U. (0.1 pc) resolution, from the 1.7-meter 
Antarctic Submillimeter Telescope and Remote Observa- 
tory (AST/RO). The large maps provide a global view 
of the entire cloud, enhancing the analysis of existing 
observations, such as those of Liseau et al. (1999) with 
ISO and large scale submillimeter continuum mapping 
observations such as Motte et al. (1998) and Johnstone 
et al. (2000). Regions of special interest are mapped at 
3000 A.U. angular resolution in CO, CS, HCO+ and their 
rare isotopomers from the 10-meter Heinrich Hertz Tele- 
scope of the Submillimeter Telescope Observatory (HHT) 
in Arizona. The distribution of all of these species, and 
their relevance to the structure of the p Oph cloud is 
discussed. In § 4, the cloud's photodissociative prop- 
erties are explored by combining CO, HCO + and [C I] 
data with a thorough far-infrared spectroscopic survey of 
the p Oph region using the Infrared Space Observatory 
(ISO) (Liseau et al. 1999). We use simple one- and two- 
component models to compute molecular column densi- 
ties and physical conditions from the sample of observed 
diagnostic atomic and molecular spectral lines. Finally 
we examine the future role of infrared H2 emission as 
a useful probe of the p Oph cloud's PDR structure and 
excitation. 

2. PRESENTATION OF OBSERVATIONS 

Several large-scale studies of the p Oph cloud have been 
performed, from millimeter-wave spectral line maps of 
C 18 0(1 -> 0) pioneered by Wilking & Lada (1983), to in- 
frared studies of the embedded stellar cluster such as Bar- 
sony et al. (1997), to (sub)millimeter continuum studies, 
most recently Motte, Andre & Neri (1998), Wilson et 
al. (1999), and Johnstone et al. (2000). Liseau et al. 
(1999) used the Long Wavelength Spectrometer (LWS) 
on ISO to obtain continuum fluxes and integrated intensi- 
ties for the far-infrared (FIR) atomic fine-structure lines 
of [C II] (158 pm) and [O I] (63 and 145 pro) toward sev- 
eral positions of the p Oph cloud and PDR. To enhance 
the interpretation of these varied studies, we mapped 
regions to maximize overlap with published data. For 
example, we chose several east-west strip positions from 
Liseau et al. (1999) for spectral line mapping. We have 
also mapped over much of the p Oph A cloud core, and 
toward denser portions of the stellar cluster. 

The observations were obtained using the AST/RO 
telescope at the South Pole, the Heinrich Hertz Telescope 
(HHT) on Mt. Graham (Arizona), the Caltech Submil- 
limeter Observatory (CSO) on Mauna Kea (Hawaii), and 
with published data from the Infrared Space Observatory 
(ISO) (Liseau et al. 1999). These data are summarized 
in Tableland described in detail below. Figure 2] pro- 
vides an overview of our submillimeter spectral line maps 
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TABLE 1 

Summary table of submillimeter observations. 



Position^' 


Line 


Telescope 


Map Size 
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FWHM 
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WL1 


CO(4 - 3) 


AST/RO 


40x40 


180 




CO(7 - 6) 


AST/RO 


3x3 


60 




[C I](l-0) 


AST/RO 


40x40 


210 




[C I](2- 1) 


AST/RO 


3x3 


60 


WL16 


CO(2 — 1) 


CSO 


5x5 


3 1 




C ls O(2 - 1) 


HHT 


5x5 


34 




C ls O(3 - 2) 


HHT 


5x5 
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HCO+(3 - 2) 


HHT 


6x5 
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H 13 CO+(3 - 2) 


HHT 


1.5x1.5 


28 




HCO+(4-3) 


HHT 


5x5 


21 




CS (5 - 4) 


HHT 


5x5 


34 




CO(7 - 6) 


AST/RO 


3x3 


60 




[C I] (2-1) 


AST/RO 


3x3 


60 




[OI](0- 1) < b > 


ISO 




62 




[Ol](l-2) < b ' 


ISO 




83 




[C H](f -|) (b) 


ISO 




68 


ISO-EW5 


C ls O(3-2) 


HHT 


5x5 


23 




HCO+(3 - 2) 


HHT 


8x5 


28 




HCO+(4-3) 


HHT 


3x6 


28 



( a ) Equatorial coordinates of central positions in epoch J2000: 
"WL1"= (16 h 27 m 01 s .5, -24°26'41"7), WL16=(16 h 27 m 02M, 
-24°37'25".5), ISO EW4=(16 h 26 m 21 s .92, -24°25'40".4), 
ISO EW5=(16 h 26 m 35M, -24°25'40".9) 
( b >Data taken from Table 2 of Liseau et al. 1999 



in the overall context of the p Oph cloud. 

Large 40' x40' maps of the p Oph dark cloud were made 
in the 492.1607 GHz (609 pm) fine-structure neutral car- 
bon line and in the 461.0878 GHz (650 pm) CO (4 -> 3) 
line with AST/RO in October 1998 (Figure 0). Our ref- 
erence off-source position (offset +9000" in RA and - 
3000" in DEC) for these observations was tested to be 
free of measurable CO line emission at these frequen- 
cies. The maps were centered near the source WL1 
(Wilking & Lada 1983) and oversampled (1' grid in RA 
and DEC). Pointed observations were performed in the 
CO(7 -> 6) and 809.3 GHz (370 pm) ( 3 P 2 -> 3 P 1 ) to- 
ward the p Oph A cloud core, the ISO EW4 position of 
Liseau et al. (1999) and the 460/490 GHz map center. 
All AST /RO observations were taken in position switch- 
ing mode and the total integration time (on and off) per 
position was 1 minute for CO (4 — * 3), 3 minutes for [C I] 
J = 1 - and 20 minutes for [C I] J = 2 — 1. The [C I] 
observations were taken with a 492 GHz quasi-optical 
receiver (Engargiola, Zmuidzinas & Lo 1994, Zmuidzi- 
nas & LeDuc 1992) in 1998 and the new 810 Ghz, 4- 
beam heterodyne array receiver, PoleSTAR in 2003. The 
CO (4 — > 3) observations were performed in 1998 with a 
461 GHz SIS waveguide receiver (Walker, Kooi and Ja- 
cobs 2001). The beam size at the time of the 492, 461, 
and 810 GHz observations was 3.5', 3' and 1', respec- 
tively. The line strengths were calibrated using hot and 
cold loads, coupled with sky measurements at the refer- 
ence position (see Stark et al. 2001). A velocity resolu- 



tion of 0.4 km s _1 was achieved for these spectra using 
a 2048 channel acousto-optical spectrometer (Schieder, 
Tolls & Winnewisser 1989). The main beam efficiency 
and system temperature for the 400 GHz receivers was « 
0.7 and 1000-2000 K, respectively. PoleSTAR's coupling 
efficiency, 0.5 (T sys ~ 20,000 K) was lower owing to very 
coarse optical alignment at the time of (first light) ob- 
servations. An instrumental velocity offset of 3 km sec -1 
at 460 & 490 GHz was removed by comparison with ob- 
served CO and [C I] lines toward infrared source WL16 
from the HHT and the Caltcch Submillimeter Observa- 
tory (CSO). 

We used the HHT to obtain overlapping CO iso- 
topomer maps (C 18 and 13 CO) for two of the east-west 
ISO strip positions of Liseau et al. (1999). Roughly 
5'x5' maps were made in C 18 0(2 -> 1), C 18 0(3 2), 
13 CO(2 -> 1), CO(2 -> 1), CS (5 — 4), HCO+(3 -» 2) 
and HCO + (4 — > 3) toward the east-west strip position 
EW4 of Liseau et al. (1999), and a limited subset of ob- 
servations toward EW5 and toward infrared source WL16 
and nearby Elias 29 (Figures I3I4|) . The maps were taken 
in March 1999, May 1999, and April 2001 using the spec- 
tral line on-the-fly mapping capabilities at the HHT with 
facility 230 GHz and 345 GHz SIS waveguide receivers. A 
2048 channel, 250 KHz resolution acousto-optical spec- 
trometer was used as a backend for both receivers, as 
were three filterbanks of 1 MHz, 250 KHz, and 62.5 KHz 
resolution. We adopted the same (clean) reference posi- 
tion (9000" ,-3000") as used at AST/RO for the higher- 
frequency observations. The half-power beamsize was 
-30" for the 230 GHz band observations and -20" for 
the 345 GHz band. The adopted beam efficiencies for 
an extended source were 0.8 for the 230 GHz receiver 
and 0.7 for the 345 GHz receiver, interpolated from mea- 
surements of point sources and the full moon. System 
temperatures with the 230 GHz receiver were typically 
400-600 K, and 500-1500 K with the 345 GHz receiver. 
Calibration uncertainty for the AST/RO, HHT and CSO 
data is estimated to be —20% below 500 GHz and —30% 
at 810 GHz. 
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[C I] and CO J=4-3 (AST7RQ 




Fig. 1. — Contextual overview of the main p Oph cloud. The grayscale image 870 (im continuum map from Johnstone et al. (2000), 
atop which the locations of infrared sources, cloud cores, and the extent of our spectral line submillimeter maps are annotated. A region 
of infrared H2 emission indicative of a PDR front is labeled by a quadrilateral. 




Fig. 2. — 40' X 40' map centered near infrared source WL1 in CO(4 — ► 3) (a) and [C I]( 3 Pi — » 3 Po) (b). For CO(4 — > 3) contours are in 
steps of 9ct=7 K km s _1 starting at 20 K km s" 1 . For [C I]( 3 Pi — > 3 Po) contours are in steps of 4<r=2 K km s _1 starting at 6 K km s — 1 . 
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Fig. 4.— - 5' X 5' maps toward position EW4. in CO(2 -> 1) (left) and 13 CO(2 -> 1) (right). For CO(2 -> 1), contours are in steps of 
15o"=5.5 K km s — 1 starting at 5.5 K km s _1 . For 13 CO(2 — > 1), contours are in steps of 9cr=3.9 K km s — 1 starting at 6.5 K km s _1 . For 
C ls O(2 — > 1) contours are in steps of 6cr=1.6 K km s _1 starting at 4 K km s _1 . For CS(5 — > 4) contours are in steps of 1.5o"=2 K km s — 1 
starting at 4 K km s" 1 . Contours for HCO+(3 -> 2) and HCO+(4 -» 3) are in steps of 2.5 K km s 1 starting at 4 K km s 1 . 
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3. MORPHOLOGY OF THE p OPH CLOUD 

3.1. Line Excitation and Radiative Transfer 

Emergent spectral line intensities from molecular 
clouds are a product of the combined effects of molec- 
ular abundances, excitation, and physical and kinematic 
cloud structure. Explicit treatment of the level popula- 
tions and transfer of line radiation is necessary to disen- 
tangle these components and obtain a meaningful picture 
of the structure of molecular clouds. To simplest approx- 
imation, we can assume that the level populations can be 
described by a Boltzmann distribution (i.e., local ther- 
modynamic equilibrium, or LTE), and determine total 
column densities in a simplified and standard way (cf. 
Keene et al. 1987): 



Nn 



1 ( Ej 

9.J \k-Lex 



Q 



T R dV 



(1) 



where J Tr dV is the main beam integrated intensity of 
the spectral line, Q is the partition function for the en- 
ergy levels accessible under these interstellar conditions, 
gj is the statistical weight of level J, and A JV is the 
Einstein coefficient for a transition from rotational state 
J — > J'. In the case of CO and its isotopomers, the 
populations of the lowest rotational levels, with nomi- 
nal n cr i t = 6 x 10 2 cm -3 in the J = 1 — > transition 
to n cr i t = 3.3 x 10 4 cm~ 3 in the J = 4 — > 3 transi- 
tion, should be thermalizcd to good approximation for 
an adopted density of n — 10 4 cm" 3 , consistent with 
ISO observations of [C II] and [O I] line emission (Liseau 
et al. 1999). 

Because the excitation of submillimeter and THz-wave 
lines of CO and related atoms and molecules may be sub- 
thermal, a more explicit calculation of the level popula- 
tions is necessary. Assuming detailed balance and steady 
state, the equations of statistical equilibrium are solved 
iteratively with the level populations and corresponding 
column densities, using the LTE level populations as ini- 
tial conditions. Using the observed linewidths, emergent 
line intensities and optical depths are computed using an 
escape probability formalism for radiative transfer (de 
Jong, Boland, & Dalgarno 1980). The competition be- 
tween radiative and collisional processes in depopulating 
a particular energy level typically involves a large num- 
ber of possible collisional transitions, the effects of ra- 
diative trapping of optically thick lines, and the rates 
of absorption and stimulated emission. We adopt the 
collisional cross sections from Schinke (1985), Monteiro 
(1985), Schroder et al. (1991), Jaquet et al. (1992), and 
Flower (1990) for H 2 interactions with CO, HCO+, C°, 
C + , and 0°, respectively, and suitably scaled for the ap- 
propriate kinetic temperature. The H 2 ortho-para ratio 
is assumed to be thermalized at the gas temperature. 
Interactions of 0° and C + with electrons are computed 
using rates from Pequignot(1990) and Blum & Pradhan 
(1992), respectively. 

A sample calculation highlights the differences in ex- 
citation between LTE and NLTE (escape probability) 
at moderate densities for various PDR- probing species. 
The results are plotted in Figure A plane-parallel 
slab of N(H 2 )= 2.5 x 10 22 cm" 2 and n = 10 4 cm~ 3 
radiates in lines of CO, C 18 0, [C I], [C II], and [O I]. 
2 km s" 1 linewidths are adopted for the transfer of ra- 
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Fig. 5. — The ratio of escape probability level populations to 
LTE populations as a function of energy above the ground state. 
Conditions: n = 10 4 cm -3 and = 25K; column densities scaled 
from N(H2)=2.5 X 10 22 cm -2 with "standard" interstellar abun- 
dances (see text). 



diation through the slab. C + and O are subject to 
T k = 50 K, the other species feel T k = 25 K. The level 
populations are compared to the predictions of LTE for 
the same total column densities as a function of energy 
above the ground state. The adopted abundances for this 
illustrative example, relative to H 2 , are [CO] = 2 x 10~ 4 . 
[C 18 0] = 4 x 10~ 7 , [C°] = 2x 10~ 5 , [C+] = 4 x 10 5 ^ 
and [O ] = 4 x 10~ 4 . 

The submillimeter fine structure lines of atomic car- 
bon are exceptionally represented by thermal excitation, 
as is [C II]. The very lowest rotational levels of C 18 
are thermalized, but the level populations rapidly de- 
part from LTE for J>5. Radiative trapping in optically 
thick low-J CO transitions maintains a thermal popula- 
tion distribution up to J<8, after which it too becomes 
subthcrmally excited. It is useful to point out that in 
cold clouds with no hotter gaseous components, the low- 
est rotational levels of CO are the only levels that are 
populated, thusly LTE can often provide a reasonable 
approximation. Our calculation of C 18 line intensities 
compare favorably with the LVG models of Mauersberger 
et al. (1992), which is expected since LVG and escape 
probability models are functionally similar. The high 
critical densities of HCO + and O lead to subthermal 
excitation in all levels above the ground state. Care- 
ful accounting of the excitation and level populations of 
these species is especially critical to the correct interpre- 
tation of their spectra. 

3.2. Morphology of Large Scale CO and [C I] Emission 

With a simplified but effective mechanism for comput- 
ing the NLTE level populations and emergent line inten- 
sities of interstellar atoms and molecules, we are prepared 
to explore the structure of the p Oph cloud. 

In Figure we present the AST/RO large scale maps 
in CO(4 -> 3) and [C I]( 3 Pi -> 3 P ), centered near in- 
frared source WL1. The peaks of the two maps differ 
by roughly one 3' beam, or 0.1 pc at the distance of the 
p Oph molecular cloud. The CO(4 — ► 3) line emission 
is very optically thick (r ~ 10 2 ) where C 18 is easily 
observed (Wilking & Lada 1983). We therefore expect 
the peak line intensity to trace the kinetic gas temper- 
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Fig. 6.— Overlay of the AST/RO [C I]( 3 Pi -> 3 P ) map 
contours atop the Abergel et al. (1996) ISOCAM image at 7 & 15 
pm. The [C I] emission is extended and generally follows the edges 
of warm dust emission expected for the illuminated surface of the 
molecular cloud. Crosses indicate cloud cores identified in FigureHl 
The thick oval highlights a region of intense H2 and PAH emission 
indicative of a H2 photodissociation front, from Boulanger et al. 
(1999). 



TABLE 2 

Thermalization of C° level populations 



Position 


Tmb (CO) 


T(C°) 




(K)( a ) 


(K)< b ) 


p Oph A 


45 


41 


WL1 


20 


22 


WL16 


24 


2V 



( b ) Derived from Gaussian fits to the peak intensity of the CO (4 
3) and (7 — > 6) lines 

( b ) Derived from both 809 and 492 GHz [C I] lines 



ature at the r ~ 1 surface of the cloud's "CO photo- 
sphere". Sample 12 CO spectra across the cloud are plot- 
ted in Figure [7| Where the CO line profile is Gaussian, 
the corrected main beam antenna temperature is gen- 
erally 10-15K although it peaks at 35K at offset (-20', 
0'), consistent with surface heating from HD 147889. In 
practice, however, this analysis is often complicated by 
prominent self-absorption in the CO line profiles toward 
the cloud cores. We fit two component (simultaneous 
emission & absorption) Gaussian functions to the 12 CO 
line profiles and recover reliable peak brightness temper- 
atures. Since we measure absorption line strengths as a 
byproduct of this analysis, we also perform a (curve of 
growth) analysis of the absorption lines to yield maps 
of the foreground cold absorbing gas. Over much of the 
40' x 40' map, the excitation of the cold gas is not well 
constrained, but certainly must be lower than the excita- 
tion temperature of the line-emitting gas. Over the range 
of 10 < T < 30K, the derived column densities vary by 
less than 30% if we assume that the cold gas is repre- 
sented by T ex =15K. The gas is certainly not colder than 
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Fig. 7.— A) AST/RO spectra of optically-thick CO(4 -* 3) 
across the p Oph molecular cloud, resampled to a 4' grid for read- 
ability. The velocity width for each spectrum is 9 km s — , and the 
corrected main beam temperature scale goes from -3 to 30K. Al- 
though the line profiles are severely self-absorbed over the majority 
of the main cloud, Gaussian fits to the line profiles yield estimates 
of the kinetic temperature and B) the structure of foreground ab- 
sorbing CO gas, derived from the degree of self-absorption. As in 
Figures fJandEl offsets are relative to the map center near infrared 
source WL1. 



10K, even in the unexcited eastern portion of the cloud, 
as the weaker CO (4 — ► 3) lines emitted there are not 
self-absorbed, as they would be if the excitation of fore- 
ground gas were markedly colder than the line-emitting 
gas. 

The CO emission peak, with a peak brightness temper- 
ature of ~45K is spatially coincident with the filamentary 
western boundary of the cloud. This strongly implies the 
foreground heating of the CO boundary due to UV pho- 
toillumination and the subsequent competition for those 
photons by H2 molecules and dust grains. The warm 
CO is localized, suggesting that the rest of the molecular 
cloud's Earth-facing "CO photosphere" is cold and not 
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subject to the significant radiation field from HD 147889 
or the other B-stars attending the dark cloud. 

In contrast to 12 CO, the overall structure of the 
[C I]( 3 -Pi — > 3 Po) integrated intensity map closely resem- 
bles the C 18 0(1 — ► 0) integrated intensity map toward 
p Oph taken by Wilking and Lada (1983) and matches 
the line intensities and profiles of Frerking et al (1989). 
The [C I] map morphology generally does not match the 
CO (4 — > 3) map, due to the large differences in line opti- 
cal depth. The addition of the 810 GHz CO (7 -> 6) and 
[C I]( 3 P2 — > 3 -Pi) lines confirms that the [C I] level popu- 
lations are thermalized at the CO excitation temperature 
(TableEl 25-45K in the west, 20-25K in the map center, 
and 10-15K to the east). The [C I] emission peak has a 
line center optical depth of < 1, similar to C 18 0(1 — > 0) 
f t!3.3.1fl . Indeed, the 492 GHz line optical depth actually 
reaches a maximum of r = 1 at the second [C I] peak near 
p Oph F due to the lower kinetic temperatures expected 
there and the correspondingly greater populations in the 
J = and J = 1 levels. 

In Figure El we overplot the same AST/RO [C I] con- 
tours atop a gray-scale image of the p Oph cloud taken in 
7 & 15 /im dust emission, from Abergel et al. (1996). The 
warm dust probed by the mid-IR highlights the heated 
"edges" of the molecular cloud, especially the filament- 
like structure at the western edge of the cloud. This 
western filament is coincident with copious rotational 
and ro-vibrational H2 line emission and PAH emission 
indicative of an H2 photodissociation front (Boulanger 
et al. 1999). The [C I] emission is clearly extended, and 
generally follows the edges of the warm IR dust emis- 
sion, as expected for a tracer of the warm surfaces of 
molecular clouds. [C I] emission peaks in two regions; 
both correspond to dense condensations of dust and gas 
seen in (sub) millimeter dust continuum intensity maps 
(Motte et al. 1998, Johnstone et al. 2000). The eastern 
intensity peak is well correlated with the dense molecular 
ridge associated with the p Oph E and F condensations 
(Motte et al. 1998) and matches a similar distribution in 
C 18 line emission (Wilking & Lada 1983). The western 
[C I] peak lies just south-east of the H 2 dissociation front 
as observed from Earth. The C 18 peak of p Oph A lies 
even further southeast of the [C I] peak, reminiscent of 
the C + /C°/CO arrangement of plane-parallel PDR mod- 
els. This suggests that some of the western PDR front is 
observed edge-on or at least obliquely. We will return to 
the structure of the western PDR front in H4.3I 

3.2.1. Relation of [C I] to [C II] and CO Line Emission 

Kamegai et al. (2003) also reports large scale, albeit 
undersampled, observations of the p Oph cloud in [C I] 
emission, and it is of interest to compare their maps 
with the present work. The previous study interprets 
the eastern [C I] line emission peak ("Peak II") as be- 
ing potentially due to the chemical youth of the region, 
because the atomic carbon and C 18 integrated inten- 
sity maxima do not follow the standard C + /C/CO ar- 
rangement expected in plane-parallel PDRs. Indeed, an 
atomic-carbon-rich region could be due to the ongoing 
formation of a dark cloud where the ionic, atomic, and 
molecular forms of carbon have not yet equilibrated. Al- 
though this prospect is exciting, our study of the p Oph 
cloud in both species does not strongly support this in- 
terpretation for the following reasons: 



1) Emission from atomic carbon does in fact trace the 
CO column density and dust continuum relatively faith- 
fully. The eastern "carbon peak" in our maps follows a 
northwest-to-southeast ridge that is also reproduced in 
C 18 J=1^0 emission in the Wilking & Lada (1983) 
map and also realized by our C 18 J=2— >1 and HCO + 
J=3^2 maps around Elias 29, part of the atomic carbon 
"ridge" (Figure 0}. Furthermore, the p Oph dust con- 
tinuum cores "E" and "F" pass through the same ridge 
(Motte et al. 1998), including sources such as Elias 29 
and a portion of the stellar cluster (see Figures ^ an d El 
for an overlay of our [C I] map and the ISO infrared map 
with dust continuum cores labeled). 

There is some discrepancy however; the C 18 ridge 
continues to the northwest, through the p Oph contin- 
uum peak "C" , whereas the atomic carbon emission turns 
to the north instead. Kamegai et al. (2003) appear to use 
such a cut through the cloud to demonstrate the differen- 
tiation of atomic carbon and CO. However, they neglect 
the fact that HD 147889 lies beyond the p Oph cloud. 
It is easy to imagine that the slightly reversed [C I] /CO 
arrangement on the sky is entirely due to projection ef- 
fects along this particular cut. Most of the [C I] emission 
follows the 13 CO or C 18 distribution accurately. 

2) A problem of timescales exists. The hypothesized 
"carbon-rich peak" in the youthful process of forming 
CO from C + and C° has already been the site of active 
star formation for ~ 1 Myr, a typical age of young stars 
in the p Oph cluster. Given the high gas densities sug- 
gested by the abundant HCO+ (3 -» 2) emission cospa- 
tial with the [C I] and CO ridge (Figure 21, the chemical 
timescales for equilibrating the CO abundance are likely 
to be appreciably shorter (< 10 6 yr) than Kamegai et al. 
(2003) assume. Their values stem from the time depen- 
dent models of Lee et al. (1996), which adopt a maximum 
gas density of only 10 4 cm -3 , over two orders of magni- 
tude less than that needed to thermalize the J = 3 — > 2 
transition of HCO + . One caveat with this interpretation 
lies in the possibility that the HCO + and atomic carbon 
trace different gas along the line of sight. A carbon map 
of higher angular resolution will help clarify the density 
structure toward this complex region. 

3.3. Structure of the p Oph A Cloud Core 

We have marked the location of our smaller scale maps 
taken with the HHT in Figure^ The centers of two such 
maps coincide with ISO-LWS pointings from Liseau et al. 
(1999) in [C II] and [O I] in an east-west strip adjacent 
to the p Oph A cloud core. A third mapped region is 
centered on infrared source WL16 (Wilking & Lada 1983) 
and includes Elias 29. 

3.3.1. CO Line Opacity and Excitation 

Figures 0] depicts the HHT spectral line maps toward 
ISO-EW4 which skirt the p Oph A cloud core and bi- 
sect a hot dust filament seen in the ISOCAM images of 
Abergel et al. (1996). 12 CO 2 1 and 4 -> 3 line 
profiles through the region are badly contaminated by 
strong self- absorption (Figure [SJ, indicating high optical 
depth and cold foreground material. A plot of r through 
the 13 CO(2 — * 1) line is shown in Figure El where a 
standard abundance ratio of 7.1 has been adopted for 

nTc^^Qi • Although this ratio is higher in the outermost 



10 



KULESA et al. 



PDR layers from poor C 18 self-shielding and subse- 
quent selective dissociation (Black & van Dishoeck 1988, 
Plume et al. 1999), when summed over the column of 
the entire cloud, this estimate is likely to be reasonable. 
We see conclusively that the 13 CO(2 — > 1) emission is 
optically thick with a line center optical depth of w 3; 
furthermore, the opacity in the (2 — > 1) transition of 
C 18 is non-negligible, about 0.5. Even modest optical 
depths lead to radiative trapping which easily thermal- 
izes the J<3 levels for densities prevalent in dense clouds 
(>10 3 cm -3 ). The 13 CO 2 — > 1 line profiles are Gaussian 
except toward the p Oph A cloud core, toward the north- 
east corner of the maps (Figure|HJ). Kinetic temperatures 
derived throughout the region are 20-35K, peaking at the 
position of the p Oph A2 core (Motte et al. 1998) and 
along the warm dust filament that marks the western 
foreground edge of the cloud. 
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Fig. 8.— HHT spectra of optically-thick CO and 13 CO in 
the (2 — » 1) transition toward the ISO EW4 position; integrated 
intensity maps may be found in Figure [I] Line profiles have been 
resampled to a 50" grid for readability. As in Figure |1] offsets are 
relative to the EW4 map center of Liseau et al. (1999). For the CO 
map, the velocity width for is 25 km s — 1 , and the corrected main 
beam temperature scale goes from -4 to 40K. For 13 CO, the velocity 
width is 15 km s — 1 per spectrum, and the corrected brightness 
temperature scale goes from -4 to 30K. 






Fig. 9. — The bottom panels portray a column density correc- 
tion for line opacity, C(t) = t/(1 — e~ T ) over observed CO and 
HCO + line profiles, derived from the ratioing of spectra of different 
isotopomers presented in the upper panels. Assumed abundances 

N( 13 CO)__-, , N(HCO+) 
arC N(C^O)-'- L and N(H«CO + ) : 



=70. 



3.3.2. Abundance and Excitation of HCO + and CS 

While examining Figure it is important to point out 
that no two tracers of molecular material appear alike in 
our maps, even those which should trace similar condi- 
tions. In the case of CO and its rare isotopomers, the 
differences result from the effects of radiative transfer 
and line optical depth, probing different depths into the 
cloud. However, HCO + and CS line emission should be 
probing gas of high density, yet the emission line maps 
share nothing in common. Nor do they sample the same 
gas seen in C 18 line emission. 

For example, the emission in the C 18 and HCO + 
transitions strongly peak toward the NE corner of our 
map, coincident with the p Oph A cloud core and proto- 
stellar source VLA1623 (Andre et al. 1993 and references 
therein). The C 18 emission in particular is roughly 
comparable to the (sub) millimeter continuum emission 
by Motte, Andre & Neri (1998), Wilson et al. (1999), 
and Johnstone et al. (2000). The C 18 line emission 
differs principally in that the p Oph A core dominates, 
and the more filamentary structures seen in the contin- 
uum do not contrast as strongly. For example, VLA1623 
is easily separated from the p Oph A core itself in the 
IRAM and SCUBA maps, but is seen merely as the end 
of an elongated ridge of C 18 line emission stemming 
from p Oph A and not distinguishable from it. 

Unlike CO and its isotopomers, HCO + emission (Fig- 
ures I10f) is not prominent throughout the map; it is 
instead largely confined to a ring-like structure that in- 
cludes the C 18 core of p Oph A at its eastern edge. 
It qualitatively reproduces the 1 mm submillimeter con- 
tinuum structure of Motte et al. (1998) (figure 2a) and 
Johnstone et al. (2000). However, the ringed, filamen- 
tary structure is more pronounced in HCO + than it is 
in submillimeter continuum emission or C 18 line emis- 
sion. It is notable that the center of the "ring" shows 
no sign of HCO + self-absorption, as demonstrated by 
the line profile maps of Figure ITUI and H 13 CO + shows 
a similar drop in intensity at the center of the "ring" . 
Two possibilities explain the heightened HCO + intensity 
around the p Oph core: [1] a heightened HCO + abun- 
dance due to enhanced cosmic ray flux, electron density 
and/or shock chemistry, [2] a normal HCO + abundance, 
but in gas that is either exceptionally warm and sub- 
thermally excited, or very dense {n ~ 10 6 cm -3 ) and 
thermalized near the CO excitation temperature. 

The column density and abundance of HCO + provides 
the best constraints on these possibilities, and can be 
determined by observations of H 13 CO + . For example, 
a line opacity analysis for the HCO + peak (EW4 offset 
+70", +70") that corresponds to VLA1623 is shown in 

Figure H assuming J|ffi^co+) = 70 - The HCO+(3 -> 2) 
emission is optically thick with r ~ 25 at line cen- 
ter. Assuming that [HCO+/H 13 CO+]=70, N(HCO+)= 
5 x 10 14 cm" 2 . To relate this value to that of CO and 
H2, the column density of C 18 is calculated by mea- 
surement of its 3 — > 2 and 2 — > 1 integrated intensi- 
ties. The derived gas temperature is 25K, consistent 
with the observed peak antenna temperature of 13 CO, 
and the column density of cold C 18 is computed to be 
1.6+0.3 xlO 16 cm" 2 , yielding N(H 2 )= 4+0.9 xlO 22 cm~ 2 
if [C 18 0]=4 x 10~ 7 (Wilson & Rood 1994) and a strongly 
enhanced HCO + abundance of 1.2 x 10~ 8 . The same 
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Fig. 10.— Extended 3' x 6' 
are spaced 3 K km s _1 apart, 
The velocity width is 10 km s~ 



map from the HHT toward the p Oph A cloud core (offsets from ISO EW4) in HCO+(4 — > 3). Contours 
starting at 6 K km s _1 , respectively. The line profiles have been resampled to a 40" grid for readability, 
per spectrum, and the corrected brightness temperature scale ranges from -4 to 15K. 
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Fig. 11. — Physical conditions derived from the integrated in- 
tensities of HCO+ 3 -» 2 and 4 3 in the escape- probability 
formalism, for the south side of the HCO + "ring" (EW4 offset 
10", 60"). Four tracks represent different abundances of HCO + , 
and the length of the tracks represents the span of parameter space 
that fit the available observations. An abundance of < 3 X 10 — 9 
best matches the column density of H 13 CO+ if the isotopomer ratio 
is 70. The excitation temperature of CO is plotted as a horizontal 
line, flanked by la uncertainties as dotted lines. 



situation holds for the peak of HCO + emission at the 
north-west corner of the HCO + "ring" (see Figure 1101 
offset 0", 180"). 

The high abundance of HCO + is localized. To illus- 
trate this point, the south side of the HCO + "ring", at 
ISOEW4 offset (10", 60") can be seen in FigureH Based 
on the detection of optically-thin H 13 CO + (3 — » 2) emis- 
sion of intensity 0.5 K km s _1 and an estimated H2 col- 
umn density of 2.5 x 10 22 cm~ 2 from the 2 — ► 1 and 
3^2 lines of C ls O, the derived abundance of HCO + is 
3 x 10 -9 , somewhat enhanced but not unusual. Further- 
more, the intense submillimcter continuum peak around 
VLA 1623 is all but absent in HCO+ and CO emis- 
sion, suggesting significant depletion, with abundances 
of HCO + not even reaching 10 -10 times that of H2. The 
parameter space of kinetic temperature and density is 



plotted in Figure El The excitation of HCO + is consis- 
tent with CO if it is subthermal at a density slightly less 
than 10 6 cm~ 3 . If HCO + stems from a colder environ- 
ment than probed by CO line emission, higher densities 
are required. 

The CS molecule is often adopted as a tracer of high 
density environments like HCO + , yet CS 5^4 line emis- 
sion in the p Oph region looks very different from HCO + . 
The p Oph A core is nearly indistinguishable from the 
high column density plateau in which it sits. The abun- 
dance of CS at the south-west corner of the filament, at 
the excitation conditions of CO and HCO + (T& = 25K) 
is only 7 x 10~ 10 relative to H2 at a density of 10 7 cm -3 . 
The derived abundances are only weakly dependent on 
excitation temperature between 20 and 50K, and no CS 
enhancement is seen in any of the regions where HCO + 
is enhanced. On the contrary, the CS abundance seems 
lowest where C 18 is concentrated. The CS distribu- 
tion peaks instead at the p Oph A2 core (Motte et al. 
1998). The abundance, given a kinetic temperature of 
35K from the C 18 measurements, is higher; 2 x 10 -9 at 
n = 10 7 cm" 3 to 3 x 10~ 8 at n = 10 5 cm" 3 . The CS line 
profile has a red shoulder, indicating that it contributes 
to a (dense) outflow in the region; indeed, the line profiles 
of 13 CO and C 18 also indicate the presence of higher 
velocity material. The enhanced CS abundance toward 
p Oph A2 is most readily explained by the liberation of 
sulphur-bearing species from dust grains to the ISM via 
heating and shocks in the outflow (e.g. Mitchell 1984). 

The differences in the CS and HCO + maps suggest that 
different chemical mechanisms are responsible for their 
relative abundance distributions. If shocks best explain 
the variation in CS, is an enhanced ionization rate a plau- 
sible explanation for the behavior of HCO + ? To explore 
this possibility, we consider a simplified set of reactions 
which govern the abundance of HCO + , namely formation 
via the ionization of H2 to (followed by H3" + CO 
— > HCO + + H2), and destruction by recombination with 
free electrons and reactions with gaseous water. The lat- 
ter is negligible due to the low abundance of water in the 
p Oph A core as measured by the Submillimeter Wave 
Astronomy Satellite (SWAS) (Snell et al. 2000; Ashby 
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et al. 2000). We adopt an electron fraction of 5 x 10 8 , 
estimated coarsely from [e] w (C/^e n {^2)) 1 ^ 2 (Langer 
1985) and increased by a factor of 3 to approximately 
account for (a low abundance of) metal species (Lee et 
al. 1996). The ionization rate is estimated by combin- 
ing a normal cosmic ray ionization rate (3 x 10~ 17 s _1 ) 
with the local X-ray flux derived from the brightness and 
projected distance from all point sources seen in public 
ACIS/Chandra images of p Oph A (Hamaguchi et al. 
2003). Three sources with Lx ~ 10 30 erg s _1 are found 
within 3' (the magnetic B star p Oph SI 2.6' away, and 
two sources within 1 , one seemingly at the position of 
VLA 1623). DoAr 21 (L x ~ 10 33 erg s' 1 ) lies at a 
separation of 6.5' and also contributes to the total X- 
ray flux (Imanishi et al. 2002). Following Maloney et 
al. (1996), the impact of these sources upon the ion- 
ization rate of (molecular) hydrogen is estimated to be 
(2 — 6) x 10~ 17 s' 1 , comparable to the contribution of 
cosmic rays. At a density of nn = 10 6 cm' 3 , this en- 
hanced ionization rate yields HCO + abundances of 10~ 8 , 
matching the peaks observed toward p Oph A. 

These chemical gradients make the p Oph A core an 
interesting chemical laboratory for the study of recent 
star formation and subsequent feedback upon molecular 
clouds that warrants further study. 

4. PDR MODELING OF THE p OPH CLOUD 

4.1. The Vicinity of the p Oph A Core 

The map that coincides with ISO pointing "EW4" 
(Liseau et al. 1999) slices through the southern-most 
portion of the hot dust filament seen in the Abergel et al. 
(1996) ISOCAM image, at the edge of a region loosely 
interpreted as the photoilluminated edge of the p Oph 
cloud. It also lies immediately adjacent to the p Oph A 
cloud core. Owing to the large amount of high-quality 
data toward EW4, we will initially focus our analysis to- 
ward this region. 

4.1.1. Column Densities 

The range of telescopes used to obtain these data re- 
sult in half-power beam sizes that vary widely between 
emission lines. To compare the derived column densi- 
ties in a consistent way, we convolve our HHT spectra 
to the resolution of AST/RO (3.5') and ISO (1.1'). An 
indirect comparison between AST/RO and ISO obser- 
vations can then be established via the convolved HHT 
maps. The convolved spectra, with #fwhm = 3.5', are 
plotted in FigurclT2lalong with the EW4 position spectra 
of CO(4 -> 3) and [C I] from AST/RO. 

The ratio of integrated intensities of C 18 in its (3 — > 
2) and (2 — > 1) transitions yields T ex = 22K for those 
levels under LTE. An escape probability calculation us- 
ing the integrated intensities of observed C 18 and 13 CO 
lines yields the column densities listed in Table [3] for 
Tkin = 30 K and n > 10 4 cm' 3 , implied by the antenna 
temperature of 13 CO 2 — * 1, the homogenous PDR mod- 
eling of far-infrared lines toward ISO-EW4 by Liseau et 
al. (1999), and the observed [O I] and HCO+ lines in 
the beam, both of which have critical densities above 
10 5 cm -3 . The resulting H 2 column density is estimated 
to be 2.1 x 10 22 cm' 2 , which corresponds to A v = 26 
using standard A V /E(B — V) = R = 3.1 dust grains, or 
A v — 18 using R = 4.5, likely to be more appropriate 
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Fig. 12.— Composite ISO-EW4 spectra from the HHT (con- 
volved to the AST/RO resolution of 3.5') plotted on same velocity 
axis as the CO and [C I] lines from AST/RO. 
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Fig. 13.— Convolved ISO-EW4 spectra at 810 GHz of CO 
(J = 7 -» 6) and [C I] 3 P 2 - 3 Pi from the PoleSTAR array at 
AST/RO. 
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Resulting column densities toward ISO-EW4. 
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to dense molecular environments like p Oph (Carrasco, 
Strom, & Strom 1973). The derived H2 column den- 
sity from 13 CO is 40% higher than that derived from 
C 18 0; however its higher optical depth, unmeasured ex- 
citation temperature and non-Gaussian line profile sug- 
gests that the C 18 observations are likely to be more 
reliable. Furthermore, the escape-probability-calculated 
optical depths of the C 18 lines are consistent with the 
previous LTE analysis that contributed Figure [5] 

The total column densities for C°, C + and 0° can only 
be calculated with knowledge of their excitation. As only 
one line is reliably observed for each of these species, we 
will now appeal to PDR models to provide the radiation 
field and excitation as a function of depth into a cloud, 
and will attempt to match all available observations. One 
can then elucidate the types of regions from which these 
three species must radiate, thereby disentangling some 
of the physical structure of the p Oph PDR. 

Liseau et al. 1999 have already performed an analysis 
of the [O I] and [C II] emission in their ISO-LWS spec- 
tra using homogeneous PDR models, and we will extend 
their results to include the CO, HCO + and [C I] results 
presented here. We will also discuss two models of sim- 
ple inhomogeneous clouds, which are offered to provide 
a better fit to all available observations. 

4.1.2. Homogeneous PDR 's 

The archetypal PDR model invokes a plane parallel 
slab of constant density, illuminated on one side by a 
source of ultraviolet photons scaled by the local inter- 
stellar radiation field, Iyy (also Go=Iuv x 2.23), and on 
both sides by the cosmic microwave background. The ba- 
sic chemical network (including PAHs) and heating and 
cooling terms (i.e. thermal structure) follows the treat- 
ment of Kaufman al. (1999). The PDR is gridded into 
zones, each with a given radiation field, density, and tem- 
perature. The treatment of the photodissociation in each 
zone is then characterized by the methods described in 
van Dishoeck & Black (1988) and Black & van Dishoeck 
(1987). The rate of H2 dissociation includes the effects 
of UV line overlap as described by Draine & Bertoldi 
(1996). Given the excitation temperature derived from 
the 63 pm [O I] transition, and based on the results of 
Kaufman et al. (specifically Figure 1 of that paper), we 
estimate a radiation field of Itjv ~ 30 for position ISO- 
EW4. We adopt a mean density of 2 x 10 4 cm -3 , which 
fits the the 63 pm [O I] and 158 pm [C II] line emission 
as per Liseau et al. (1999), with a slab of total col- 
umn density N(Ha) = 2 x 10 cm" 2 . Our models predict 
the correct line intensity for [C I] 3 Pi — > 3 Po emission, 
2.3 x 10 -6 erg s" 1 cm" 2 sr , though the Kaufman et 
al. (1999) models overestimate the [C I] flux by 50%. 
Both models similarly fit the 158 pm [C II] and 63 pm 
[O I] line intensity within the observational uncertain- 
ties. As described by Liseau et al. (1999), the 145 pm 
[O I] emission is underproduced by an order of magni- 
tude in the models. Furthermore, the observed intensity 
of HCO + (3 — > 2) and (4 — > 3), with critical densities 
of ~ 2 x 10 6 and ~ 1 x 10 7 cm" 3 , respectively, is not 
reproducible at n ~ 10 4 cm" 3 unless the excitation is 
exceptionally hot (T > 300K) and subthermal. Such 
conditions are not indicated by the relative intensities of 
the HCO + (3 — > 2) and (4 — ► 3) lines, which are more 
consistent with a thermalized population at T« 20K, 



similar to CO. Higher densities are simply necessary to 
create sufficient submillimctcr HCO + emission. A natu- 
ral way to account for the discrepancy in excitation is to 
adopt a heterogeneous model for the physical structure 
of the region. 

4.1.3. Clumpy PDR's 

The simplest heterogeneous models for PDR's involve 
the consistent composition and weighting of two homoge- 
neous PDR models. That is, we consider a plane-parallel 
slab of gas containing dense condensations (clumps) of 
gas randomly embedded in a more diffuse envelope, and 
illuminated by an ultraviolet radiation field. The dense 
clumps individually receive an attenuated radiation field 
based on the amount of intervening diffuse material pho- 
tons must travel to reach the clump. Scattering in the 
interclump medium is neglected. We consider two cases: 
one in which the clump:envelope density contrast is 30:1 
(10 5 and 3 x 10 3 cm" 3 ), and a second in which the con- 
trast is 1000:1 (10 6 and 10 3 cm" 3 ). The clumps are as- 
sumed to be spherical, 0.1 parsecs in diameter, with a 
volume filling factor of 0.1 and 0.2, respectively, so as to 
keep the mean density comparable to the homogeneous 
models, n ~ 10 4 cm" 3 . The clump sizes are consistent 
with the assumed sizes for clumps in other multicompo- 
nent PDR models (e.g. Meixner & Tielens 1993; Spaans 
1996) and correspond to the angular size of major C 18 
and HCO + condensations in our spectral line maps. 

The two inhomogenous PDR models give similar qual- 
itative results, which are summarized in Figure PHI The 
large majority of the [C I] 3 Pi — * 3 Po emission comes 
from the diffuse envelopes, as has been suggested by 
Spaans (1996) and MT93, among others. The envelope 
exceeds the low critical density and excitation energy 
needed to produce significant [C I] line emission, and 
as it dominates the filling factor of the region, it also 
dominates the production of [C I] 3 Pi — » 3 p). This is 
especially true of the high-density-contrast model. The 
excitation of the [C I] is expected to be similar to that 
of 13 CO and C 18 0, which is borne out by the similar 
distribution of [C I] and C 18 in p Oph, as well as the 
similarities in line opacity and profile. 

The somewhat higher excitation temperatures and 
critical densities of [C II] emission make its point of ori- 
gin slightly more ambiguous. Almost 2/3 of the [C II] 
emission still arises in the more diffuse envelopes, but 
the remaining line emission stems from the warm surface 
layers of the denser clumps, particularly those closest to 
the PDR front. 

Low-J 12 CO line emission is highly saturated, and 
hence comes from the r ~ 1 surface layer of the molec- 
ular cloud "photosphere" at the frequency of each CO 
line, and is equally excited in the diffuse envelope and 
the surfaces of the clumps that are near the PDR "sur- 
face" . Since the emergent line intensity is weighted by 
filling factor, most low-J line emission comes from the dif- 
fuse envelope. As one moves into the sub-millimeter and 
far infrared bands, the less-saturated high-J lines arise 
most prominently from the warm, dense clump surfaces. 
Line emission from the isotopically rare species, espe- 
cially C 18 0, is also dominated by the envelopes, but to a 
lesser extent than 12 CO, as the reduced self-shielding of 
the rarer species in the UV-dominated transition regions 
leads to isotope-selective dissociation and fractionation. 
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TABLE 4 

Homogeneous PDR model results for ISO-EW4 



Diagnostic 


Observed 


Modeled 


Reference 


[C II] 158 fim 


2.4 X 10" 4 < a ) 


2.5 x 10" 4 < a > 


Kaufman et al. (1999) 


[O I] 63 Aim/ [C II] 158 urn 


0.43 


0.4 


This paper( b ) 






0.9 


Kaufman et al. (1999) 


[C I] 3 Pi - 3 P 


2.3 x 1CT 6 < a > 


3.5 x 10" 6 M 


Kaufman et al. (1999) 






2.4 x 10" 6 


This paper( b > 


CO(2 - 1) / [C II] 3 P 3/2 - 3 Pi/ 2 


~ 50 


100 


Kaufman et al. (1999) 


HCO+(4 -* 3) 


5.9 K-km s" 1 


0.36 K-km s" 1 


This paper( b ) 



( a ) ln erg s 1 cm 2 sr 1 

( b ) Cloud model follows that of van Dishoeck & Black (1988), see text 



TABLE 5 

INHOMOGENEOUS PDR, DENSITY CONTRAST = 30: 3 X 10 3 cm~ 3 

& 10 5 cm" 3 



Diagnostic 


Observed 


Modeled 


[C II] 158 Atm 


2.4 x 10" 4 W 


2.0 X 10~ 4 ( a > 


[O I] 63 fim/ [C II] 158 /im 


0.43 


0.36 


[C I] 3 Pj _> 3 P 


2.3 x 10~ 6 ( a > 


2.1 x 10~ 6 ( a > 


CO(2 -> 1) 


1.2 x 10" 6 W 


1.0 x 10" 6 W 


HCO+(4 -+ 3) 


5.9 K km s" 1 


1.6 K km s" 1 



' a 'ln erg s 1 cm 2 sr 1 



TABLE 6 

INHOMOGENEOUS PDR, DENSITY CONTRAST = 1000: 10 3 cm -3 h 
10 6 cm" 3 



Diagnostic 


Observed 


Modeled 


[C II] 3 P 3/2 - 3 P 1/2 


2.4 X 10 -4 (a) 


1.8 x 10~ 4 (a) 


[O I] 63 Ami/ [C II] 158 /an 


0.43 


1.0 


[C I] 3 Pi -» 3 P 


2.3 x 10" 6 W 


2.8 x 10" 6 W 


CO(2 -» 1) 


1.2 x 10~ 6 < a > 


1.0 x 10~ 6 ( a > 


HCO+(4 -> 3) 


5.9 K km s" 1 


5.2 K km s" 1 



( a )ln erg s 1 cm 2 sr 1 



The lower line optical depth in these transitions probes 
the depth of the entire cloud, not just the surface layers 
seen in low-J 12 CO. 

Far-infrared [O I] 3 P — > 3 Pi and submillimeter HCO+ 
line emission is only efficiently produced at higher den- 
sities, so it doesn't significantly contribute to the line 
emission in the diffuse envelope, but instead stems from 
the dense clumps. The HCO + (4 — » 3) and (3 — > 2) emis- 
sion is readily excited at low temperatures, and there- 
fore stems from the cold, shielded, molecular clump in- 
teriors as well as the UV-illuminated PDR surfaces of 
clumps. The relatively high 228K excitation energy of 
the [O I] 3 Po — ► 3 Pi transition limits the vast major- 
ity of the intense [O I] emission to the warm surfaces of 
the clumps. However, if atomic oxygen is not entirely 
bound up in molecules like CO, H2O, or depleted onto 
dust grains in the shielded cores of dense clouds, it can 
emit at detectable levels even at relatively cold tempera- 
tures (T fc = 10-20 K) at high densities (n ~ 10 6 cm" 3 ). 



Indeed, it now seems certain that O2 is not a significant 
repository of elemental oxygen in molecular clouds (cf. 
Goldsmith et al. 2000, and references therein); however, 
the nature of oxygen depletion in molecular clouds is still 
debated enthusiastically. 

A narrow range of density contrast is required to simul- 
taneously account for the modest strength of the 63 /im 
[O I] and submillimeter HCO + emission. At the higher 
densities (n ~ 10 6 cm -3 ) needed to stimulate significant 
HCO + (4 — > 3) emission, [O I] emission can become very 
efficient and rapidly overwhelms the 158 fjxa [C II] inten- 
sity for normal interstellar conditions, contrary to ISO- 
LWS observations of p Oph where [C II] emission is domi- 
nant (Liseau et al. 1999; Kaufman et al. 1999). However, 
clumps with nonuniform (power law) density, n(r), are 
theoretically expected and naturally satisfy both excita- 
tion conditions. As such, submillimeter HCO + emission 
would then originate primarily from the dense cores of 
the clumps (n > 10 6 cm -3 ), whereas the [O I] emission 
would stem from the warmer "surfaces" of the clumps 
(n w 10 5 cm" 3 ). 




Fig. 14. — Illustration of a clumpy PDR and the physical region 
each emission line probes. 



4.2. Other ISO positions 

Our large scale maps of [C I] overlap with all 13 point- 
ings of ISO-LWS from Liseau et al. (1999). We can 
therefore perform limited (homogenous) PDR modeling 
toward a 1-dimensional east-west cross section of the p 
Oph cloud. In Figure IT51 we plot the cross-sectional in- 
tensities of 158^m [C II], 63pui [O I] and 609/mi [C I] 
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line emission. As described in Liseau et al. (1999), 
the [C II] line emission decreases as projected distance 
from HD 147889 increases, tracing the incident radiation 
field upon the p Oph PDR. With a high critical density 
and excitation energy, the [O I] flux traces both density 
and radiation field intensity. [C I] line emission however 
shows a different spatial relation, instead following the 
cold gas column density as traced by C 18 (for compar- 
ison, see Wilking & Lada 1983); the two maxima in the 
[C I] distribution occur as the cross-sectional slice skirts 
the edges of the p Oph A and p Oph B2 cloud cores. 

to 

O p-i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i-q 




RA Offset (arcmin) 

Fig. 15.— Intensities of 158/xm [C II], 63^m [O I] and 609Aim 
[C I] line emission along the ISO east-west cross-section of the p 
Oph cloud 

As in § 14.1.21 we use the ISO-LWS integrated intensities 
of 158 pm [C II] and 63 pm [O I] line emission to estimate 
Itjv and n. The thickness of the PDR model is decreased 
with increasing projected distance from HD 147889 in 
order to correctly recover the correct [C I] line intensi- 
ties. Model results for four representative points along 
the east- west cross section are presented in Table [7] 

4.3. Directly Probing the Western Hi PDR Front 

Mapping across the western PDR front in the infrared 
emission lines of H2 can clarify the cloud geometry at its 
"true surface". Boulanger et al. (1999) presents ISO H2 
emission line spectra from the western PDR front of the 
p Oph cloud in the 0-0 S(0) through S(7), and 1-0 S(l) 
transitions. The H2 emission is observed to be coincident 
with both PAH emission in the thermal IR and the fila- 
mentary dust structure seen with ISOCAM at 7 and 15 
pm. The non-thermal population of high H2 rotational 
levels coupled with photoilluminated PAH emission sug- 
gests that UV fluorescence, rather than thermal shocks, 
are likely to be responsible for the H2 excitation. Fluo- 
rescing H 2 emission directly probes the incident UV field, 
PDR density, and geometry of the H I/H 2 transition re- 
gion at A v < 1 with high angular resolution. 

Spectroscopy of fluorescent H2 emission can also elu- 



cidate the location of the H2 dissociation region rela- 
tive to the rest of the molecular cloud. For example, 
the H2 1-0 S(0,1) and 1-0 Q(2,3) line pairs arise from 
the same upper states. As these are both intrinsically 
weak electric-quadrupole transitions, the effects of self- 
absorption can be safely neglected. Thus the relative in- 
tensities of these two lines are functions only of their rel- 
ative spontaneous radiative transition probabilities and 
photon energies. For example: 

Iq(3) = ^Q(3) Ag(i) ^ 

h(l) ^Q(3) 

for the quadrupole H2 A-values computed most recently 
by Wolniewicz, Simbotin & Dalgarno (1998). 

In Table line ratios are tabulated for visual extinc- 
tion (A v ) from to 40. Such measurements, combined 
with the intensity of the emission, would not only con- 
strain the physical conditions in which H2 is dissociated, 
but also would clarify whether the UV-illuminated sur- 
face that is fluorescing in H2 is on the front or back side 
of the p Oph cloud. 

For example, at Ttjv = 50 and n ~ 10 4 cm" 3 , the 
expected H2 1-0 S(l) line intensity is modest, about 
3.5 x 10" 6 erg s" 1 cm" 2 sr" 1 . If the observed 1-0 S(l) 
emission (from Boulanger et al. 1999) originates from the 
back side of the cloud, as does the rest of the [C II] emis- 
sion, the near-IR H2 emission would have to pass through 
a mean column density loosely equivalent to A v ~ 25, as- 
suming a mean reference value of the gas/dust extinction 
ratio of 1.6 x 10 21 cm" 2 mag" 1 measured in diffuse in- 
terstellar clouds. This translates to a mean extinction of 
~ 3 magnitudes at 2.12 pm, interpolating the Rieke & 
Lebofsky (1985) infrared extinction curve. This would 
extinguish the H2 1-0 S(l) line flux by a factor of ~ 17, 
to 2.1 x 10~ 7 erg s" 1 cm" 2 sr" 1 . Noting this faintness, 
and pending photometric calibration of the H2 1-0 S(l) 
images of Boulanger et al. (1999), we predict that the H2 
emission instead comes from a region of molecular ma- 
terial on the front side of the cloud that is not shielded 
from HD 147889 and therefore is subject to photodisso- 
ciation. A schematic view of the p Oph PDR might then 
look like Figure IT^l 

Sensitive, high resolution near-infrared long slit spec- 
troscopy of H2 in PDRs provides a natural complement to 
the kinds of long wavelength PDR diagnostics descibed in 
this paper, and are worth renewed effort with numerous 
far-infrared and submillimeter facilities soon to become 
available (such as Herschel & SOFIA). 

5. SUMMARY 

In this paper, we have attempted to exploit the 
complementary information afforded by large-scale and 
high-resolution submillimeter spectroscopic observations 
made with small and large antennas. 

We present large scale maps (40' x 40'), obtained with 
AST/RO, in CO (4 -> 3) and [C I]( 3 Pi -> 3 P ) centered 
near the infrared source WL1 in the p Ophiuchi molecu- 
lar cloud. [C I]( 3 Pi — » 3 Po) emission is extended on large 
scales, resembles the C 18 column density map presented 
by Wilking and Lada (1983), and generally follows the 
edges of warm dust emission expected for the illuminated 
surface of the molecular cloud. The optically-thick CO 
(4 — > 3) emission traces the CO "surface" of molecular 
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TABLE 7 

PDR Modeling of a p Oph cross-section; n = 10 4 cm -3 



RA Offset^) 


Ajv 


TpDR 


N(H 2 ) 


63/im [O I]/ 


609Aim [C I] 


(from WL1) 




(K) 


(cm -2 ) 


158/mn [C II] 


(erg s — 1 cm -2 sr — 1 ) 


-18'=EW1 


90 


200 


6 x 10 21 


0.4 


1 X 10" 6 


-6'=EW5 


30 


100 


1 x 10 22 


0.5 


2 x 10" 6 


3'=EW8 


20 


100 


1 x 10 22 


0.2 


1.5 x 10" 6 


18'=EW13 


10 


50 


4 x 10 21 


0.1 


8 x 10" 7 


( a ) "EW" positions correspond to Liseau et al. (1999) 



TABLE 8 

H2 V = 1 — > LINE RATIOS AS A FUNCTION OF INCREASING 
EXTINCTION 



A v 


lQ(3)A3(l) 


T (a) 





0.703 


3.5 x 10" 6 


5 


0.785 


2.0 x 10" 6 


10 


0.875 


1.1 x 10" 6 


20 


1.09 


3.7 x 10~ 7 


30 


1.36 


1.2 x 10" 7 


10 


1.70 


3.9 x 10" 8 



'^Extinguished intensity of 1-0 S(l), in erg s 1 cm 2 sr 1 




: 120 pc 
t 

/^^^ observer 

Fig. 16. — Schematic diagram of the p Oph dark cloud, with 
PDR surfaces highlighted in black. The concave back side of the 
cloud, per Liseau et al. (1999), bears the bulk of the UV radiation 
field from HD 147889. However, if the extinction toward the H2 
emission (Boulanger et al. 1999) is low, the "filament" PDR on 
the western edge of the dense molecular cloud likely represents a 
foreground PDR region that is not shadowed from HD 147889 by 
the rest of the cloud. The depicted dark cloud is also surrounded 
by a diffuse "halo" of gas and dust. 



clouds and is strongest (hottest) at the western bound- 
ary of the molecular cloud, adjacent to exciting B2V star 
HD 147889. Adjacent ro-vibrational H2 emission, cospa- 
tial with a western filamentary structure in the hot dust 
continuum at 7 & 15 pm and in PAH emission, provides 
evidence for a new foreground PDR. This is in contrast to 
the rest of the cloud, which appears to be illuminated on 



the far side by HD 147889. Additional H2 spectroscopy 
is needed to verify these suppositions. 

We also present smaller scale maps (~ 5' x 5'), taken 
with the HHT, in C 18 0(2 -> 1), C 18 0(3 -> 2), 13 CO(2 -> 
1), CO(2 -> 1), CS(5 -> 4), HCO+(3 -» 2) and 
HCO + (4 — > 3) toward a single region (EW4) in the 
AST/RO map, chosen to overlap with [C II] and [O I] 
observations made by Liseau et al. (1999) with ISO- 
LWS. The C 18 structure of the p Oph A core is directly 
comparable to recent submillimetcr continuum emission, 
however HCO + line emission is concentrated in a ring- 
like structure. Modeling of the HCO + ring in comparison 
to CO suggests that the condensations are due to both a 
higher HCO + abundance and a dominant dense gas com- 
ponent with n ~ 10 6 cm~ 3 . A slightly enhanced ioniza- 
tion rate contributes to the elevated HCO + abundance 
in this region. Although it also stems from dense gas, 
the emission from CS follows a totally different morphol- 
ogy. The p Oph A cloud core contributes to a plateau of 
CS emission, but the actual CS peak is coincident with 
the relatively minor p Oph A2 core and outflow instead. 
Shock and outflow chemistry is suspected to contribute 
to a higher CS abundance in this region. 

We use the emission lines presented in this paper to 
infer column densities of all optically thin species to- 
ward ISO-EW4, and three other positions examined by 
ISO-LWS. We use the indirect determinations of H2 col- 
umn density from C 18 and 13 CO with PDR models 
to constrain the physical conditions in which [O I] and 
[C II] line emission is observed (Liseau et al. 1999). 
We model the physical structure of the p Oph PDR in 
the traditional homogeneous, plane-parallel formalism, 
and obtain self-consistent fits to most available observa- 
tions. However, the luminous and clumpy submillime- 
ter HCO + and infrared [O I] emission directly implies 
dense substructure. We adopt a simple two-component 
PDR model, featuring a diffuse gaseous envelope in which 
dense clumps reside, to account for the HCO + and [O I] 
emission. In this picture of the physical structure of the 
PDR, we are able to attribute emission from different 
species to different physical locations within the p Oph 
PDR. We find that [C II] and [C I] emission predomi- 
nantly arises from the lower density envelopes. We find 
little evidence for non-equilibrium chemistry providing a 
"carbon-rich" region as per Kamegai et al. (2003), find- 
ing instead that the atomic carbon distribution follows 
the C 18 ridge through dust continuum peaks E and F. 
We find that like [C I], low- J CO emission stems predom- 
inantly from the envelopes, and high- J emission should 
stem from the denser clumps. Submillimeter HCO + and 
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infrared [O I] emission stems from the clumps and clump 
surfaces, respectively, and indicate clump surface tem- 
peratures of 50 — 200 K, an ultraviolet radiation field 
with I uv - 10 - 90 (G ~ 10 2 ), densities of 10 5 to 10 6 
cm~ 3 , and interior temperatures of < 20K. 

This study illustrates the use of far-infrared, sub- 
millimeter, and infrared spectroscopy for dissecting the 
complicated structure of clouds and star forming re- 
gions. Yet, numerous new advances in detector tech- 
nology will soon enhance our understanding more dra- 
matically, putting the evolution of clouds and stars in 
a more global, Galactic context. The advent of large- 
format infrared detector arrays will make broad and nar- 
row band observations of hot dust and gas possible over 
entire clouds, as well as providing photometry for de- 
tailed extinction maps. Although integrated heterodyne 
detector arrays at submillimeter wavelengths are in their 
infancy, rapid progress is expected to revolutionize the 
radio universe as arrays have already done for the in- 
frared. Finally, new airborne and space-based platforms 



such as the Spitzer Space Telescope, Herschel, SOFIA, 
and SAFIR will open up the far infrared, opening new 
windows into the chemical history of carbon, nitrogen, 
oxygen and water in the Galaxy. 
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